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Abstract. 

Some open questions exist with fluctuation-induced forces between extended 
dipoles. Conventional intuition derives from large-separation perturbative approxi¬ 
mations to dispersion force theory. Here we present a full non-perturbative theory. In 
addition we discuss how one can take into account finite dipole size corrections. It is 
of fundamental value to investigate the limits of validity of the perturbative dispersion 
force theory. 
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For many years there has been strong interest, both theoretically and 
experimentally, in fluctuation induced forces m- Until now most work on dispersion 
forces using either the quantum electrodynamical or the semi-classical formalism have 
relied upon a point dipolar description. The relationship between these two theories 
has been discussed by Ninham, Parsegan, and Weiss [6j. Previous work that indeed did 
incorporate the finite size of an extended dipole was based on perturbative expansions 
(truncated to the leading term), suggesting that the van der Waals energy contributed 
with a short-range attractive binding energy [7j. Our non-perturbative approach has 
been investigated for atom-atom dispersion forces in water [8j and for extended dipole 
self-energies. [9j The effects of a non-perturbative theory on self-energies (which can 
give 50% changes or more compared to perturbative theory) ought in principle to be 
detectable experimentally, although it is a challenge to measure them directly. Solvation 
energies of atoms in a dielectric medium, i.e. changes of the self-energy in a vacuum 
compared to in a medium, can be measured unmu. Use of an incorrect theory may 
lead to problems in interpretations of experimental solvation energies. Conceptually 
important additional effects of the non-perturbative theory will be demonstrated in 
the present work for dipole-dipole dispersion interactions in vacuum. Perturbative 
expansions in terms of multiple dipole-field interactions are often restricted to leading- 
order approximations where each dipole scatters the field only once. This is a valid 
approach if the coupling is very weak and/or the separation between the dipoles is 
large enough. The aim of this work is to demonstrate how the keeping of the full 
non-perturbative, multi-scattering theory and taking the finite spread of the dipole into 
account, alters the non-retarded van der Waals force at contact distances. The inclusion 
of finite dipole size in a non-perturbative theory opens up for the possibility of having 
van der Waals repulsion when two dipoles come very close to each other. 

In the following we first recapitulate the theory of van der Waals energy between 
point like dipoles. Then we introduce a theory that accounts for finite spread of dipoles 
in the van der Waals force between two dipoles in their ground states. The van der Waals 
potential is created by a sum of symmetric and antisymmetric modes. These modes can 
be related to the non-retarded first order dispersion forces. [12j Repulsion in the van der 
Waals force may occur as some modes stop contributing when the dipoles come close 
together. Within this theory avoiding perturbative series expansions, we find substantial 
corrections to the dipole self-energy. [9] We also consider the resonance interaction when 
one dipole is in its ground state and the other is in an excited state. At resonance the 
excitation switches back and forth between the two dipoles. The interaction can then 
be separated into three branches, the j = x-, y- } and z-branches, where the name of 
each branch denotes in what direction the oscillating dipoles are pointing. We let the z- 
direction be defined as the line joining the two dipoles. We focus on deriving analytical 
asymptotes that can be used for all separations within the non-retarded regime. The 
correct way to evaluate resonance interaction in the retarded limit is still debated [121113] . 
We present results in the non-retarded limit, ft is reassuring that in this limit different 
groups obtain the same results. [T2UT5] As we will demonstrate when two dipoles come 
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close together the results remain finite in contrast to the textbook results (which is a 
l/p 3 -dependence for resonance interaction and l/p 6 -dependence for ground state van 
der Waals potentials where p is the interdipole separation). 

Mitchell et ah published in 1972 a paper [5] demonstrating short-range repulsion 
between two point-like dipoles when they come very close together. Consider two 
identical dipoles with polarizability given by 

a{iuj) = a 0 /{l + uj 2 /ujl), (1) 

with Uq being the ionization frequency of the dipole. The allowed coupled frequencies 
in the limit of small separations are given by Eq. (27) in Mitchell et al. [5]: 

w = w 0 \/l± a 0 /p 3 ,ui = u 0 y/l± 2a 0 /p 3 , (2) 

p being the distance between the dipoles. The first two roots each occur twice. An 
important point is that only real frequencies contribute to the interaction energy. If all 
roots are real, i.e. 2ao/p 3 < 1, then the interaction energy is 

U{p) = 

x l 2 ^ 1 + «o/p 3 + 2a/1 - a 0 /p 3 (3) 

+ \/l + 2«o/p 3 + a/ 1 — 2«o/p 3 — 61. 

For large distances this expression goes over to the London interaction energy. 
The key point here as demonstrated by Mitchell et al. [5] was that, for decreasing 
separations, first 2 a 0 /p 3 > 1 an d then a 0 /p 3 > 1, several modes become imaginary and 
do not contribute to the (real) energy. The interaction energy becomes in this limit 

U(p) = |2a/1 + a 0 /p 3 + \/l + 2a 0 /p 3 — 6 j 

~^{(2 + V 2 )^/ p 3 / 2 - 6 }. 

The freezing out of modes at distances of the order of the dimensions of the dipoles 
leads to a repulsive energy. Of course the dipole approximation itself breaks down in 
this limit. 

This problem is partially relaxed by including spread out size effects of the extended 
dipole. A difficulty still occurs, however, as can be seen from the fact that the logarithm 
in the usual expression for the interaction energy between two identical molecules, for 
example, turns imaginary when x > 1 for some of the Matsubara frequencies. One has 
to sum over terms of the form ln(l — x). The problem will be addressed in the following 
discussion. 

From the equations of motion for the excited system it is straightforward to 
derive the zero temperature Green’s function for two identical dipoles in vacuum 
HHEO. The resonance condition [12] can be obtained from the following relation: 
1 — a(tu) 2 T(p\uj) 2 = 0, where 1 is the identity tensor, p is the dipole-dipole separation 
distance, T is the susceptibility tensor and a(u>) the polarizability (at frequency co) of the 
dipole. In the case of two identical dipoles excited in the jth branch the above resonance 
condition can be separated into one antisymmetric and one symmetric part. Here we 
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refer to the quantum mechanical state function. In the symmetric case the state function 
is the sum of the function where dipole 1 is excited and dipole 2 not and the function 
where dipole 2 is excited and dipole 1 not. In the antisymmetric case the state function 
is the difference between these two states. Since the excited symmetric state has a much 
shorter life time than the antisymmetric state the system can be trapped in an excited 
antisymmetric state [13j . The resonance interaction energy of this antisymmetric state 
can be evaluated by a simple expression for two identical dipoles, [12] 

OO 

Uj(p) = h J ^ ln [l + a(/0'/;,,(p if)] ■ (5) 


The corresponding van der Waals (Casimir-Polder) interaction between identical 
dipoles in their ground states is given by the following expression: 


OO 

u cp(e) = \ V J [i - aK) 2 y(p|*0 2 ] 


( 6 ) 


j=x,y,z ^ 
J —oo 


Let us pause here to discuss the interaction energy in terms of electromagnetic 
normal modes. The interaction energy at a specific separation is the sum of the zero 
point energy of all modes minus the corresponding sum at infinite separation. Some of 
the modes are attractive and some repulsive. Most of the interaction is canceled out 
but a net attraction survives. When the dipoles are brought closer together the energies 
of the attractive modes decrease while those of the repulsive modes increase. In the 
formalism associated to each mode corresponds one with minus the same energy; these 
negative-energy modes do not contribute to the interaction. When the dipoles come 
close enough together the energy of the most attractive mode approaches zero value. 
When this happens the one with minus the same energy reaches zero from the negative 
side. With further reduction of the separation the two modes move away, one along 
the positive and one along the negative imaginary axis. This may happen for several 
modes before the dipoles are in direct contact with each other. Our interpretation is 
that when this happens the modes stop contributing to the interaction. Since it is only 
attractive modes that drop out, the attraction gets weaker and the net interaction can 
become repulsive. 

With our interpretation one could perform the calculation in one of three equivalent 
ways: One is the mode-summation method. This method can only be used in rare cases 
when the distinct modes can be found. Then the energy of the modes that reach the 
zero value should be put equal to zero. This holds when the separation continues to 
decrease. Another way is to proceed as we did above. Then the integration should be 
performed not at the imaginary axis but along a vertical line an infinitesimal distance 
from and to the right of the axis. The result is real valued since the imaginary part of 
the integrand is odd and the real part an even function of f. A third way is to integrate 
along a closed path encircling the whole positive real axis. All three ways give the same 
result. 
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Perturbative theory. 

The traditional way to treat the integrands of Eqs. ([5]) and © is to make a series 
expansion of the logarithm and keep the lowest order term only, i.e. ln(l + x) ~ x. 
This is permissible for large separations when the lowest order term is small compared 
to unity. In this section, we consider this perturbative case while in the next section, 
we refrain from doing this approximation and term it the non-perturbative case. The 
second improvement we do is to take the finite size of the dipoles into account instead 
of treating the dipoles as point particles. The polarization cloud of real dipoles has a 
finite spread and we consider as an interesting case a spatial distribution of the dipole 
following a Gaussian function im We have, following the formalism developed by 
Mahanty and Ninliam j!4j . derived the Green’s function elements that account for 
retardation, background media, and finite dipole Gaussian radius (a) in an accurate 
way. We now consider the effect of the finite dipolar size on the non-retarded resonance 
interaction between two dipoles in an excited state in vacuum. In this limit (velocity of 
light —>■ oo) we can use the non-retarded Green’s function elements, 


and 


= 7=7- [v^Ferf® -2(f) e 



(7) 


_ -2 
\drp 3 


-0Ferf (£) + 2 (f) 




( 8 ) 


We begin by studying the effect of the finite size when we expand the logarithm 
and truncate to leading order in the dipole-field interaction. We find for the resonance 
interaction with the x-branch excited the following long range result in the non-retarded 

limit, 


u™(p) = -[2fi/(V?P 3 )] 


X 


, v 2-1 oo 

V^ erf (£) ~ 2 (£) e“(») f §a(i£). 


The corresponding result with the z branch excited is, 


(9) 


(7™(p) = [4fi/( v ^p=)]{ v ^erf(f)^ 

-2(f) ll + (f)1 e- ( -G /#««)• 


( 10 ) 


It is of interest to present the corresponding results at close contact. Taking t = p/a 
the p dependence for the x-branch is contained in fit ) = [y / 7rerf(t) — 2t exp(—t 2 )]/t 3 . 
The small t expansion of the error function is erf(t) = (2/ y/H) exp(— t 2 )(t+2t 3 /3+0[t 5 ]), 
and therefore /(f) = 4/3 + 0[t 2 ]. This produces the limit T xx —> —4/(3i/7ra 3 ). One 
can in the same way show that T zz —> T xx . Therefore in the limit of small dipole- 


dipole separation we find U Z R = U X R —> 3 ~^ 3 f ^a(if). As a result, unlike the 

well known long distance behavior [T3] . when two dipoles come very close together 
the resonance interaction is the same for all excitation-branches. We observe that the 
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resonance interaction at close contact depends on the radius oc a~ 3 . This should be 
compared to the corresponding result for ground state van der Waals interaction that 
depends on the radius oc a~ 6 . The analytical asymptotes are only valid for dipolar 
systems where the perturbative, single-scattering expansion exists. Finite size effects 
soften the attraction at very small separations. For isotropically excited dipole pairs the 
first expansion term in the logarithm cancels out and the leading non-retarded term is 
oc a~ 6 . 


Non-perturbative theory 

Now we move to the results of the non-perturbative theory. The complete 
expressions for the non-retarded first order dispersion potentials follow trivially when 
Eq. (J5]) is combined with Eq. (|7J) and Eq. (JEJ) . The simple expressions found for the 
finite size non-retarded resonance interaction suggest that one can solve this problem 
analytically, at least if one chooses a simple model polarizability (a(i£) = a(0)/(l + 
£ 2 /(Uq)). We then find the following expression (valid when the dipoles are close enough 
together so that retardation effects can be neglected): 


dx In 


x! Q (°) y r 


/„■ 

= fc„Re [-1 + Jl + a(0 yrjp 


1+x 2 


( 11 ) 


The non-retarded Casimir-Polder potential between two dipoles with finite radii is 


USp(p) = ! T E Re 

j=x,y,z 

+v /l -“(°>E™ ■ 




( 12 ) 


Asymptotically the interaction approaches the long range non-retarded limit with a 
p _6 -dependence. In the limit of zero separation between two dipoles we obtain, 


Ug?(p = 0 ) 


^Re 


—2 T y/l + 4cr(0)/(3\/7rn 3 ) 


+ y/l - 4a(0)/(3 1 /7ra 3 ) 


(13) 


This energy limit for two overlapping polarization clouds is made up of a bonding 
and an antibonding self-energy contribution which should be compared with the self¬ 
energy of a single dipole given later in this work. While the van der Waals repulsion 
between point like dipoles goes towards infinity the previous expression provides the 
upper limit to the van der Waals repulsion between finite sized dipoles. This dipole- 
dipole contribution is part of the “attractive” term in the Lennard-Jones potential. Our 
simple expression for the non-retarded van der Waals potential can easily be included in 
simulations that are using Lennard-Jones potentials. Only real modes contribute to the 
interaction energy [5j . It is important to note that besides the attractive term one should 
also add the repulsive term in the Lennard-Jones potential corresponding to exchange 
interactions due to wave-function overlap and higher-order multipole transitions PSHTZj. 
For the perturbative case, using the formalism of Richardson [18] we found in an earlier 
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Table 1 . The finite size van der Waals binding energy, U^p(p = 0) using the 
expression derived here, and as comparison the dispersion self-energy 0, U S , for 
noble gas atoms. The superscript ’truncated’ indicates that the result is from using a 
truncated (or perturbative) expansion of the logarithm in the integrand. All energies 
are in eV. 


Element 

t/“' (0) 

^truncated (Q) 

if" 

jj truncated 

He 

29.06 

-409.6 

71.21 

131.53 

Ne 

58.20 

-1046 

104.95 

220.28 

Ar 

8.767 

-133.7 

37.56 

62.07 

Kr 

4.366 

-87.98 

29.94 

47.50 


work |8j that the dipole-dipole interaction gives larger contribution than the dipole- 
quadrupole and quadrupole-quadrupole interactions at all separations. Multipolar 
contributions may, however, provide as much as 50% of the dispersion energy of some 
atoms, electron-rich anions in particular HU- They have been found necessary for 
reproducing experimental properties such as ion solvation energies ra, entropies and 
partial molar volumes (2UJ, surface tensions |2T| and activity coefficients [22j. Finite 
size corrections prevent the attraction from going to infinity when two dipoles come 
very close together. In the limit of small dipole-dipole separations the van der Waals 
interaction from the full expressions may turn into a constant repulsive energy when the 
polarization clouds overlap, but there are also effects for two dipoles at contact distances 
that reduce the attractive binding energy. The interaction approaches asymptotically 
the well-known attractive p _6 -dependence. Perturbative expansions have also been used 
to calculate the dispersion self-energy of a single dipole in vaccurn. [T]. A more complete 
theory gives the following dispersion self-energy [9], 


U s = ^ 


-1 + \/l + 4a(0)/(3a 3 y / 7r) 


(14) 


When series expanding and truncating this expression we can rederive the dispersion 
self-energy found by Mahanty [7]. However, the validity of a series expansion assumes 
that a(0)/a 3 is much smaller than unity which is not always the case. We provide 
in Table I the finite size van der Waals binding energy and the finite size dispersion 
self-energy for a number of atoms. The input data were taken from Refs. [9| [23l)24] . 
We give the results found when using our non-perturbative theory and as comparison 
also the approximate theory using a truncated series. The corrections to the self-energy 
of atoms are substantial and measurable. In fact, the non-perturbative self-energies 
in our examples are up to a factor of 2 different from the results from the truncated 
theory. These effects of a non-perturbative theory ought in principle to be detectable 
experimentally. They should be accounted for to avoid incorrect interpretations of 
experimental solvation energies and permeabilities across membranes. 


Conclusion 
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We conclude that finite size effects keep the interaction finite when two dipoles come 
close together. Notably the resonance interactions between two identical dipoles are the 
same, and attractive, for x-, y- , and z-excitation branches when the two atoms come 
close enough together. The Casimir-Polder interaction energy which is attractive at long 
range approaches a constant repulsive van der Waals energy at short range. At distances 
beyond a few dipole radii the different interactions approach well known results from 
perturbation theory. The non-perturbative theory is conceptually important, but it is 
also important for the binding energy of dipole pairs and self-energies. We emphasize 
that the reason why the interaction stays finite when the dipoles come close to each 
other is the spead out of the polarization cloud. 

There are a number of potential applications for a non-perturbative theory of 
intermolecular forces. Take for example the solvation energies of ions in a dielectric 
medium (i.e. changes of the self-energy in a vacuum compared to in a medium) that can 
be measured [XDJ[TT] . Latimer et al [25 J were able to fit experimental solvation energies 
to the Born equation by increasing the effective radius of the ions. Self-energy changes 
have also been shown to influence permeabilities across membranes [26] . The results 
obtained for the permeability of atoms across a membrane will be changed due to changes 
in the self-energy in vacuum. A factor of 2 difference for these self-energies compared 
to those obtained from a series-expanded theory ought therefore to be measurable with 
existing experimental equipment for solvation free energies and permeabilities. 

MB and CP acknowledge support from the Research Council of Norway (Project: 
221469). PT acknowledges support from the European Commission. SYB gratefully 
acknowledges support by the German Research Council (grant BU 1803/3-1) and the 
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